On the Full Column-Rank Condition of the 
Channel Estimation Matrix in Doubly-Selective 

MIMO-OFDM Systems 

Abstract 

Recently, this journal has published a paper which dealt with basis expansion model (BEM) based least-squares 
(LS) channel estimation in doubly-selective orthogonal frequency-division multiplexing (DS-OFDM) systems. The 
least-squares channel estimator computes the pseudo-inverse of a channel estimation matrix. For the existence of the 
pseudo-inverse, it is necessary that the channel estimation matrix be of full column rank. In this paper, we investigate 
the conditions that need to be satisfied that ensures the full column-rank condition of the channel estimation matrix. 
In particular, we derive conditions that the BEM and pilot pattern designs should satisfy to ensure that the channel 
estimation matrix is of full column rank. We explore the polynomial BEM (P-BEM), complex exponential BEM (CE- 
BEM), Slepian BEM (S-BEM) and generalized complex exponential BEM (GCE-BEM). We present one possible 
way to design the pilot patterns which satisfy the full column-rank conditions. Furthermore, the proposed method is 
extended to the case of multiple-input multiple-output (MIMO) DS-OFDM systems as well. Examples of pilot pattern 
designs are presented, that ensure the channel estimation matrix is of full column rank for a large DS-MIMO-OFDM 
system with as many as six transmitters, six receivers and 1024 subcarriers. 

Index Terms - BEM, full column rank, ICI, MIMO, OFDM and rank-nullity theorem. 
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On the Full Column-Rank Condition of the 
Channel Estimation Matrix in Doubly-Selective 

MIMO-OFDM Systems 

I. Introduction 

A basis expansion model (BEM) based [ 1 ] least squares (LS) channel estimation scheme for doubly-selective 
orthogonal frequency-division multiplexing (DS-OFDM) system was proposed by Tang et al in @. For good 
performance, the channel estimation matrix T needs to be of full column rank [2 pp. 2232]. However, the full 
column-rank condition was not addressed in [2, pp. 2232]. The purpose of this paper is to derive the conditions on 
the BEM and pilot patterns such that the V is of full column rank. Four BEMs are analyzed, namely, polynomial 
BEM (P-BEM) [2], complex exponential BEM (CE-BEM) 0, Slepian BEM (S-BEM) (3]- 0] and generalized 
complex exponential BEM (GCE-BEM) (2). 

A BEM-based channel estimation scheme for doubly-selective multiple-input multiple-output OFDM (DS-MIMO- 
OFDM) system was proposed by Song and Lim in |5). The methods of |2l and J5] are similar in spirit although Q 
did not deal with MIMO systems. Full column-rank conditions on the channel estimation matrix were not addressed 
in 0. In this paper, the full column-rank conditions is extended to doubly-selective MIMO OFDM (DS-MIMO- 
OFDM) systems as well. One possible way to design pilot patterns, that satisfy the full column-rank conditions, 
is presented in this paper. The pilot-pattern design examples are applicable to DS-MIMO-OFDM systems as well. 
To the best of the author's knowledge these pilot-pattern design examples have not been presented before in the 
literature. 

Tang et al authored a subsequent paper [7] which deals with the full column-rank condition of "P. The context 
and scope of the results presented in this paper differ with those presented in [7|. The key differences between the 
analysis in this paper and J7J can be summarized as follows: 

• The work of Q assumed a frequency-domain Kronecker delta (FDKD) pilot pattern. The analysis in this paper 
is not restricted to any particular pilot pattern. The FDKD pilot pattern is treated as a special case in this paper. 

• The work in (7) showed that the FDKD pilot pattern ensures that "P is of full column rank for CE-BEM and 
P-BEM. In this paper, we derive general conditions that a BEM should satisfy for V to be of full column 
rank when any pilot pattern is used. Four BEMs are investigated, namely, P-BEM, CE-BEM, S-BEM and 
GCE-BEM. 

• We also present a detailed approach to pilot-pattern design that enables V to be of full column rank. This 
design example is extended to DS-MIMO-OFDM systems as well. Furthermore, our analysis show that the 
design of pilot patterns can be independent of the conditions that a BEM should satisfy. 
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• This paper presents full column-rank conditions for DS-MIMO-OFDM systems. The work in |7| is limited to 
DS-OFDM systems. 

• Where possible, we also show that some of the results derived in this paper is in agreement with the ones 
derived in 0. 

Cao and Li presented some full column-rank conditions of channel estimation matrix for DS-MIMO-OFDM systems 
J6] Lemma 1]. We show that some of the results arrived at in this paper are in agreement with [6, Lemma 1]. 

The paper is organized as follows. The BEM-based OFDM system model is presented in Section II. This forms 
the basis for the full column-rank conditions discussed in Section III. Pilot-pattern design examples for DS-OFDM 
and DS-MIMO-OFDM systems are presented in Section IV, while Section V deals with miscellaneous conditions 
and consensus of results with prior art. Simulation studies are presented in Section VI and Section VII concludes 
the paper. 

A. Notation 

A boldface large and small letter mean a matrix and a vector, respectively. The quantity (x) denotes x mod 
N (N is the number of subcarriers in the OFDM system) while (x) a denotes a; mod a. (.) H , () T and * denote 
Hermitian, transpose and conjugate, respectively. The (p, q)th and pth elements of A and x are denoted by A(p, q) 
and x(p), respectively. The element-wise product of two matrices, called as Hadamard product, is denoted by 0, 
i.e., the (z, j)th element of A B is A(i, j)B(i, j). Ijv represents an NxN identity matrix. Where possible, we 
use the Matlab notation with the exception being vector and matrix indices, which start from zero. The notation 
B = e a Q means that B is a matrix whose dimension is the same as that of Q and B(i,j) = e aQ ^' j \ We will 
use A(:, m:n)(A(m:n, :)) to extract the submatrix from column (row) m to column (row) n, A(r, c) to extract a 
submatrix within A defined by the index-vector of desired rows in r, and the index-vector of desired columns in 
c. The vector x(m:n) extracts entries from m to n and x(r) extracts entries denoted by the index-vector r. The 
NxN FFT matrix F is defined as F(m, n) = e~~ « . The rank of a matrix A is denoted by r(A). The vector 
space spanned by vectors Xi,... ,xjv, is denoted by £(xi, . . . , xat). The vector space spanned by the columns of 
A is denoted by C(A). £(A)_L£(B) means that the column spaces of A and B are orthogonal to each other, i.e., 
A H B = 0. is an all-zero matrix or a vector whose dimension is understood from the context. S(x) denotes the 
Kronecker-delta function, i.e., 6(0) = 1 and 6(x) = for x ^ 0. 

II. System Model 

A. OFDM System Model 

Consider an OFDM system with N subcarriers. Without loss of generality, we consider the zeroth OFDM symbol. 
The N x 1 OFDM symbol vector at the transmitter, which consists of data and pilot symbols, is denoted by x. The 
kh channel tap at the nth time instant is denoted by hi(n). The length of the cyclic prefix (CP) is L samples 
(which is also the number of distinct resolvable multipaths of the DS channel). The samples at the receiver, r(n), 
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are corrupted by additive white Gaussian noise (AWGN), w(n). After discarding the CP, the iVxl received vector 
can be written as 

f = — LHF H x + n (1) 

VN 




where f = [f(L), . . . , f(L + N — 1)] T , n = [w(L), . . . ,w(L + N - 1)] T is the AWGN vector and H(p,g) = 
h/ p ^ q \(L + p) is the time-domain channel matrix J2). An FFT is applied to r to yield 

(2) 

where y = [y(0), . . . , y(N — 1)] T is the vector of demodulated output symbols. The input-output relation is therefore 
given by y = Hx + n- 

B. BEM Channel Model 

Let the /th-path channel vector be denoted as h; = [hi (L) , . . . , hi (L + N — 1)] T . This can be quite accurately 
characterized by a Qth order BEM as 1 2 1 

hi « Bh ; (3) 

where h; = [ftj(O), . . . , hi(Q— 1)] T is the vector of BEM parameters and B is an NxQ BEM matrix. Note that the 
N samples of the channel are characterized by Q BEM parameters and Q <C N. 

C. Pilot and Data Clusters 

The system has Np pilot clusters and Nd data clusters, equispaced and interleaved with each other as shown 

in Fig. Q] The lengths of data and pilot clusters are Lry = 2wrj + 1 and Lp = 2wp + 1, respectively. We denote 

the pilot and data vectors by the Nxl vectors p and d, respectively. The quantity p(i) denotes the pilot symbol 

mapped to the ith subcarrier and is equal to if the ith subcarrier is a data subcarrier. Likewise, d(i) denotes 

the data symbol mapped to the ith subcarrier and is equal to if the ith subcarrier is a pilot subcarrier. Note that 

x = p + d. Equation (O can be rewritten as 

1 - b 1 - h 1 
y = — FHF p + — FHF d + -^Fn- (4) 



Let Pf, and P scp denote the index of the middle subcarrier of the zeroth pilot cluster and the distance between the 
middle subcarriers of two neighboring pilot clusters, respectively. Since the pilot clusters are equispaced it follows 
that N = NpP sep . Let pi be an NpLpX 1 pilot-index vector whose entries comprise the indices of all pilot 
subcarriers in ascending order, i.e., 

pi =[{P b -w P ),...,P b ,..., (P b +w P ),..., (P b +(Np-l)P scp +w P }f- (5) 

The pilot cluster consists of two parts. A centrally located observation cluster (OC) of length = 2B C + 1 is used 
for channel estimation. It is protected on either side by a guard band of length G which prevents interference from 
the neighboring data clusters. Its length is dependent on the Doppler spread of the channel. The pilot cluster length 
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is therefore given by Lp = 0+2G. In all, V—NpO pilot subcarriers are used for channel estimation. The V pilot 



subcarriers are accessed via (2B c + l) pilot-index vectors pj , i — — B c ,. . ., B c , where p { 2 l> is the ith pilot-index 
vector. It comprises of the indices of the ith neighbour of the middle subcarrier in each pilot cluster and is defined 



f = [{P b +i) , (P b +i+P sep ), . . . ,{P b +i+(N P -l) P aep )f, \i\ < B c 



P2 



(i) 

Concatenating all p;, into a master pilot-index vector p 2 



(6) 



(-Bo)T 
P 2 



(B C )T 
,P 2 



, the observation pilot symbol 



vector used for channel estimation is the Vxl vector y = y(p2)- 



D. BEM-Based OFDM System Model 

Consider the P;,th demodulated pilot symbol that belongs to the zeroth pilot cluster. Neglecting interference from 
data clusters and AWGN, it follows from © that 

i=w P L—XN—1 

•v- y first,..., (iVp — 1) st pilot cluster 



1=0 n^O 



zeroth pilot cluster 

where Af—l/N. Rearranging and interchanging the summations, we have 



L— 1 wp 



N-X 



-j2n(P b +i)Afl 



hi (L + n) e l2mAfn - 



(8) 



1=0 i=—w P 



n=0 



zeroth pilot cluster 

Define the inverse fast Fourier transform (IFFT) matrix of h; as 



first,. ..,(./Vp — l)st pilot cluster 



H 



(0 



(l/iV)F ff h ; 



(9) 



The above equation represents the Doppler spectrum of the Zth channel tap. Equation ([8]) can then be rewritten in 
a concise form as 



^p {P b -wp) e-^ p "- w ^ A f l U ( p\N-w P ) + --- + p{P b + w P ) e -jMn-t«>p)A/i H a) (wp ) 
y {Pb) ~h> +P (P b +P scp -w P ) e-^( p »+ p ^-^ A f l H ( p(P scp -wp) + ■■■ 



(10) 



Similarly, the (P b + P scp )st demodulated pilot subcarrier of the first pilot cluster, can be written as 

L1 p (P b + P scp -Wp) e -j2-(P 6 +P S cp-^p)A/Z H (0 (TV — Wp) + • • • 

y(P b +P scp ) = J2 +p(P b +P scp +wp)er^ p ^^ A f l U^(wp) (11) 
' = ° +p (P b +2P scp -w P ) e -iMiW ser ^P)A// H (0 (Psop _ wp) + . 

Before we proceed further, we present a physical interpretation of ( TTOb and (fTTb . The quantity 
p (P b — wp) e~ : > 2 ' !T ( Pb ~ Wp ) A f l Tl F ^(N— wp) represents the response of the Zth path of the channel at Subcarrier P b 
due to the excitation by Subcarrier (P b — wp). The quantity e~i 27r ( Pb ~ Wp ) A f l is due to the delay of the Ith path 
channel. The quantity Hp\N—wp) represents the spillage/leakage of the kh path Doppler spectrum. It can also 
be viewed as the intercarrier interference (ICI) contribution of Subcarrier (P b — wp) on Subcarrier P b via the Ith 
channel. It therefore follows that ( TTOb and ( fTTb can be viewed as the sum of the frequency-domain convolutions 
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Fig. 1. An example of pilot and data clusters. 



of the Doppler spread of each time-shifted channel path with the pilot cluster. A detailed mathematical discussion 
and pictorial depiction of this interpretation was given by Scneiter in |8l . 

Observe that (TV — wp) is associated with p (p, — wp) in dTOb . while in ( fTTT > it is associated with p (Pf,+P scp —wp). 
Let the row-vector pc = [p(A + i-Psep — wp), ■ ■ ■ , p(Pb + iPscp + wp)] denote the iih pilot cluster. The pilot- 
pattern matrix, P pa t. is obtained by stacking all the pilot-cluster row vectors and Q is a pilot-indices matrix whose 
ith row comprises of the indices of the ith pilot cluster, i.e., 

P> - w P , 



pat 



D (0) 




Pc 






, Q = 


„(JVp-1) 




Pc 





P b + {N P - 1) P m 



Wp, 



Pb + w P 



Pb + {N P - 1) P scp + w P 



(12) 



For the example in Fig. Q] the matrices P pa t and Q are given as 

p(0) p(l) p(2) 



pat 



p(8) p(9) p(10) 



Q 



1 2 
8 9 10 



We are now in a position to present the expression for y(p| )■ To do so, we define the following 

e — j2irQA/J 
e -j2«QAfl q 
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R 



(i) 

I 



Qi (i : N P + % - 1, :) 



T>(0) t,(1) 



R 



(Afp-1) 
I 



F ff ((iY+( Pl (0)-P fc -*)),:) 
F*«JV+( Pl (l)- A -*)>,:) 



F H «7V + (pi(^ Pj L p - 1) - P b - i)) , :) 
[ho r hLj' 



(13) 



Observe that the rows of the matrices R> , R, , 



(0) R (l) j^(JVp-l) 



i.e., for instance we have 



R[ 0) (1,:) = R, (1) (0,:) 



are circularly shifted with respect to one another, 

(14) 
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Furthermore, we observe that 

rW (at — o„„^ tt(0 



H£ (JV - , . . . , H$> (w P ) , (P scp - w P 



W (0) h, ~ W (0) Bh, 



(15) 



v = 



P (A - w P ) e-JMft-wWl 

p (P 6 + P scp - W) P ) e -j27r(P i ,+P 8ep - W p)A/i 



p (P b + P sep - top) e -3^(P b +P s ^-w P )Afl 
p (Pfe + 2P sop - W P ) e -J'27r(P i ,+2P 8ep - W p)A/i 



(16) 

As already discussed, < fT3T > represents the Doppler spectrum of the (time-aligned) Zth channel tap. Recall the fact that 
interference from data clusters and noise is neglected for the moment and y(p^) = [y(P&), y(P& + P S cp), • • • , ] T ■ 
It follows from dTU), CD, 0, and CD that 



y(p(°))=H(pW : JP = 



h 



Considering the interference due to data clusters and AWGN, recalling that y = y(pa), we have 

y = Ph + m 

where ni is the interference due to data clusters and noise and 

o W(- B °)B ... L _iW(- s ^B 



(17) 



(18) 



flnW^B 



e L _iW( fl ''B 



4>L-1 



(19) 



is the channel estimation matrix of dimension VxLQ. In the above equation, <pi is associated with the Zth multipath 
and has the form 

r e l 



4>! 



W (- B ^B 

w^b 



(20) 



The purpose of this paper is to determine the conditions on the BEM matrix B and pilot clusters pi such that P 
is of full column rank. 

The matrix P defined in ( fT9l and the one defined in J3] (14)] are permuted versions of each other and have the 
same rank. The interchange of rows (columns) of a matrix does not change its rank. Equation ( fT8l is a permuted 
version of (14) in (2|. More specifically, y in this paper is defined as 

y = [y(P 6 -P c ), . • .,y(P 6 + (^p-l)P scp -P c ), . . . ,y(P 6 ), . . .,y(P 6 +P c ), . . . ,y(P 6 + (^ P -l)P scp + P c )] T 

(21) 

while yW in E (14)] is defined as 



[y(A-Pc), . . .,y(P b +P c ), • • . ,y(P 6 + (^p-l)P scp -Pc), • . . ,y(P 6 + (7Vp-l)P scp +P c )] 



(22) 
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Observe that (|2TT > and (l22l are permuted versions of each other. Similarly, h in this paper and h in [2j (4)] are also 
permuted versions of each other. 

III. Full Column-Rank Condition of V 

Recall that V is a V x QL matrix and the full column rank necessitates that r('P) = QL. The matrices 
0i;O<i<L — 1, are of dimension V X Q. It therefore follows from ( fl9l ) that 

r(&) = Q, 0<i<L-l, 

£(0,) J_ C(4h), i^j 

needs to be satisfied for V to be of full column rank. To satisfy the second condition in the above equation, it 
follows from (f20b that the relationship 

0f0,=O, ijtj (24) 

needs to be satisfied. We now present the condition that needs to be satisfied by a BEM to ensure that T is of 
full column rank. The dimensions of the matrices 0/ and E are V x VLp and VLp x Q, respectively. Generally, 
we have the relationship V 3> Q- Recall that <pi = #/E. To ensure that the first condition in d23l is satisfied, it 
therefore becomes necessary that a BEM satisfies 

r (E) = Q (25) 

which will be termed as the BEM condition (BEMC). The BEMC necessary condition for the matrix V to be full 
column rank means that the BEM matrix B containing the independent Q basis functions should note lose rank 
when it is projected in the Doppler (frequency) domain, where the observations are collected. 

Before we proceed further, we present the rank-nullity theorem. Consider a matrix B such that r (B) = Q and 
bo, . . . , bg_i, be the basis vectors that span the column space of B. If none of bo, . . . , bQ_i, lie in the null space 
of another matrix A, then the rank-nullity theorem [ 1 1 states that r (AB) = Q and Abo, ■ ■ ■ , AbQ_i, span the 
column space (basis vectors) of AB. 

IV. Pilot-Pattern Design 

The design of pilot patterns should be such that d23l i and d24"i i are satisfied. In this section, we outline one possible 
method of pilot-pattern design, that satisfies the conditions in ( 1231 ) and d24l ), thus ensuring the full column-rank 
condition of the channel estimation matrix. The design example is presented for both DS-OFDM and DS-MIMO- 
OFDM systems. 

A. DS-OFDM Systems 

We now design the pilot pattern matrix, P pa t, such that "P is of full column rank for a DS-OFDM system. 
Initially, we concentrate on the design of pilot patterns such that the second condition in d23l is satisfied. Towards 
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the end of the section, we concentrate on the design of pilot patterns such that the first condition in ( |23l > is also 
satisfied. We begin by examining R[ ' . From ([TZi i. it follows that the rth column of Q can be written as 

Q(:,r)=Q(:,0)+rl (26) 

where 1 is a column vector with all entries equal to Unity. It therefore follows that r(Q) = 2. Recall from (flTt that 
R^ = 0; (0 : Np — 1, :) = e~J 27I 'Q A /' P pat . Now consider what happens when all entries of P pa t are equal to 
Unity. It follows from [\3[ and d26l i that all columns of ; (O : Np — 1, :) in this case are scaled versions of each 
other, i.e., 



% (0:Np- 1, 



-j2ttQ(:,0)A/; -j2irrAfl -j2nQ{:,0)Afl 



(27) 



Define an N P x 1 vector as f« = [l, e -^Mi ; . . . ; e -j27riA/i(iVp-i)j ^ where A ^ _ (i/jV P ). The quantity 
can be considered as the ith harmonic of a Np x 1 vector whose fundamental discrete frequency is A/i. We 
observe that f<» ff f<>') = 0, i ^ j, < i,j < N P - 1. Recall that the rows of Q differ by P scp , N = N P P scp and 
Af = l/N. From the definition of Q in O, it follows that 

e -j27rQ(:,0)A/i = e -j2-K{P b -w P )Aflf(l) 

(28) 

£ / e -j27rQ(:,0)Am = £ . 

Under the assumption that all pilot symbols are equal to 1 (entries of P pat is equal to Unity), it follows from (|27| > 
and (f28ll that 

r(R< 0) ) =1, £(Bf ) =£(f(0) • 



Furthermore, recall from (fl4] > that the rows of the various matrices, R, , < i < Np — 1, are all circularly shifted 
with respect to one another. It therefore follows that 

r (Rp>) = 1, £ (R ; (0) ) =•••=£ (R^-^) = C (fW) • 
Next, we design P pat such that the following conditions are satisfied 



£ (r/' ) ... £(R, V ' 0; L-\ 



(29) 



£(R|°))±...±/:(Ri l 1 ) 

which would ensure that (|24] > will be satisfied. This ensures that the second condition in d23l l is also satisfied. The 
matrix t is the concatenation of the matrix, e -i 27r Q A /' q P pat , with itself, one on top of another. From (fT2l . ( fT3b 
and d28l ), it follows that Z (O : N P — 1, :) can be written as 

R| 0) = 0,(0 : iV F - 1 >: ) = [ e -^Q(:,o)A/ iPpat (:>0) _ _ _ ^Q^-DAfip^ { .^ Lp _ x) ] . (30 ) 
At this juncture, we recall the following: 

• It follows from ( fT3] > that the rows of the various matrices R, , < i < Np — 1, are all circularly shifted with 
respect to one another, as was depicted in ( fT4t . This is depicted in the upper part of Fig. [2] 

• If the columns of P pat are chosen from fW,0 < i < Np — 1, then the first condition in ( 1291 is eas- 
ily satisfied. We explain this by way of an example. Let P pat (:,0) = f^. It follows from ( fT3l l. 



9 



that R.[ 0) (:,0) = 0,(0 : N P - 1,0) = e -o^(Ph-w P )Afi f similarly, R, (1) (:,0) = 0,(1 : A^ P ,0) = 
e -i2 W A/ le -i2^(F 6 - W p+p„ p )A/; f (i + i) > Hence it foUows that £(r| 0) (:,0)) = £ (r, (1) (:, 0)) = £(f< i+1 )). 

(i) 

This holds good for other columns of R ; , < i < Np — 1, as well. This is depicted in the lower part of 
Fig. El 

In order to satisfy the second condition in d29l i. we need to distribute the vectors f W among the columns of 
Ppat- What should be the mapping strategy? Define £(a,b) = £({( a \{( b ' > ) and let P pat (:,r) = f(^). We observe 
that £(a)-L£(a),a ^ b. It follows from (O that 

Ro 0) = £(go,...,gL P -i) 

R| 0) = £(l + g ,....l + g Lp -i) (3D 



R^-! = £(L-l + <? 0l ...,£-l + <7ip-i) 



The above equation essentially says that the column space spanned by the rth column of R ; is 

£(r| 0) ) =£{l+g r )- 

In all, the L matrices Rj , < I < L — 1, have LpL columns. To satisfy the second condition in 
becomes necessary to select g r ,0 < r < Lp — 1, (select columns of P pa t) such that all the LpL columns of 
R| o) ,0 < 1<L-1, are orthogonal to each other. In other words, 

£{h+g rl )±£(h+g r2 ), 0<h,l 2 <L-l, 0<r u r 2 <L P -l, (l u n) ^ (l 2 ,r 2 ) ■ 

The above condition can easily be satisfied by the mapping 

g r = rL 

which implies that designing P pa t as 



(32) 



(33) 



pat 



f (0) f (i) _ _ f ((L p -l)L) 



ensures that 



C(0i 



£ 



£ R 



(Np-l) 
I 



£ (f( i ),f( i+ '', • • ■ , f {(Lp-1)L+1)\ 



(34) 



(35) 



Consequently, the design as per ( f34b ensures the full column-rank condition of "P. From the above equation, it 
follows that the column spaces of Matrices Rj ' , < I < L — 1 , are interleaved with one another, as is depicted 
in Fig. [3] 

We now concentrate on the first condition in ( f23l , i.e., r (<pi) = Q, < i < L — 1. Matrix E is associated with 
the choice of BEM. The quantities 0;, < I < L — 1, are associated with the pilot-pattern design and <pi — 0/E. 
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(<>) 



in 



Col. 



Col. 1 



Circularly shiftedversions 
of one another 



If columns of V are chosen from f 



R 



Col. Col. 1 



R 



Col. Col. 1 



span the same column 
space 



Fig. 2. Properties of Matrices R ; W , < i < N P - 1 




Part of column space spanned by Rj, ' 
Part of column space spanned by 



i(f< M >) L(f (L) )L(f (i+1) ) Z^f' 2 "*) 




Part of column space spanned by 



(0) 



Fig. 3. Column spaces of Matrices Rf ) , < i < L - 1, are interleaved with one another. 



It will be shown in the appendix and in Sec. VI, for the BEM models considered in this paper and the design of 
pilot pattern as per d34| i, none of the basis vectors of E lie in the null space of 0/. < I < L — 1, i.e., 

6ie Q ^0,q = 0,...,Q-l, I = 0, . . . , L - 1- (36) 

This condition will be called as rank-nullity BEM condition (RNC-BEM). This implies that if the BEMC and RNC- 
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BEM are satisfied, the rank-nullity theorem ensures that the first condition in ( |23l l, i.e., r (</»j) = Q, < i < L — 1, 
is also satisfied. It therefore follows that the pilot-pattern matrix designed as per d34l i ensures "P is of full column 
rank. In the appendix, we derive the RNC-BEM conditions in more detail for DS-OFDM systems. The derivations 
for DS-MIMO-OFDM systems can be pursued in a similar manner and is therefore not dealt in this paper. 

Note that the Np x 1 vector f W is the ith harmonic of a complex exponential vector with fundamental frequency 
A/i = (l/Np). The various f(0 constitute the columns of P pat and these ensure that £ (0/) = L (f W, f ( L + l ) i . . . j((L P -i)L+i)y 
In all, L P L vectors fW,i = 0, . . . , L P L - 1, are associated with the L matrices 9i, I = 0, 1. In order to 

ensure that 0/^0, = 0, i 7^ j, is satisfied, it becomes necessary that fW-^fw = 0, i 7^ j, < i, j < LpL — 1, and 
for this to be possible it becomes necessary to have 



B. FDKD Pilot Patterns in DS-OFDM Systems 

All pilot symbols in a FDKD pilot cluster, except the middle one, are equal to zero. Furthermore, the non- 
zero pilot symbols of all pilot clusters are assumed to be the same and equal to \. It then follows that P pa t(: 
, (Lp — l)/2) = \xi • - • , x] T ■ All other columns of P pa t comprise of entries equal to 0. It follows from ( fT2] >, ( fT3l 
and d2Sl> that 



This ensures that the second condition in d23l and (l24l l are satisfied. However, a key requirement for this to be 
possible is that f! m ) ff fW = 0, < m, n < L — 1. This implies that we must have 



Furthermore, note that r {61) = 1, r (0/) = 2B C + 1. If the BMC and RNC-BEM are satisfied, the first condition 
in (|23| ) is satisfied if 

2B C + 1 > Q- 

More discussion on the above equation is presented in Sec. V. 

C. DS-MIMO-OFDM Systems 

The system model and the channel estimation matrix, T, is an extension of the one derived in [2 ] and differs 
from the one derived in ^ in two aspects, namely, ;) the channel estimation scheme in [5] considered the case 
of O = 1 only, while the schemes in this paper and [2| consider O > 1, it) as per the scheme of (5] (6), (7)], a 
demodulated pilot symbol belonging to the ith pilot cluster takes into account the contributions of pilot symbols 
only in the ith pilot cluster, whereas, the schemes presented in this paper ( fTOb and [2| take into account the effect 
of all pilot clusters. 

Without loss of generality, we assume that all quantities correspond to the zeroth MIMO receiver. Furthermore, 
we assume that there are Nt transmit antennas. Let denote a vector of BEM parameters associated with 



L P L < N P ■ 



(37) 




(38) 



N P >L- 



(39) 
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the Zth path between Transmitter i and Receiver 0. Define — 



> n L-l 



as the concatenation of all 



hj f, I = 0,...,L— 1. The quantity P^Jt is defined in a way similar to that in (fT2l i and denotes the pilot-pattern 
matrix used by Transmitter i. Replacing P pa t in ( Tl3b by PLJ t , one obtains 0^. Using 0, W in Qj), one obtains 



R 



(fe,j) _ fl (0 



0, w (k:N p + k- 1), 



(i) 



R (o,i) R (l,i) 



pat' 

(JVp-l,i) 



. The channel estimation matrix associated 



between Transmitter i and Receiver is denoted as "P^ where 



W W (-B C ) B 



W iW (-B c ) B 



W iW (B c ) B 



, 



« 



W (- S c) B 

W^) B 



(40) 



The quantities uf' , R[ fe '^, 0,^ and are associated with the Zth path between Transmitter i and Receiver 0. 
The system model (fT~8T > can be rewritten as 

y = 7'h + ni (41) 

where V = V {a \ "P^" 1 ) , h = [h(°) T , . . . , h.^ NT ~^ T ] T and ni is the interference due to data clusters 
and noise. For the case of DS-MIMO-OFDM systems, d23l can be rewritten as 

r(0, (i) ) = Q, 0<1<L-1, 0<i<N T -l 



(<#) ± Cffi), (Mi)?* toy- 



Similar to (1241 . Ppa t , < % < N T - 1, should be designed such that 

9l? B e%>=o, (i,h)^UM 



(42) 



(43) 



When the BEMC is satisfied and none of the basis vectors of E lies in the null space of ) , < I < L — 1, < 
i < Nt — 1, the first condition in d42l > is satisfied. If we design P^l t such that 



pW = 

pat 



f (i_L P L) r^LpI+L) r(iLj»I/+(Xp-l)L) 
i ,i , . . . , i 



(k i) (i) 

it follows from the definitions of R£ ' and 0) that 



C 



(0«)=£( 



R 



= •••=£ ('r^' p_1 '^ > ) = £ ({(iLpL+l) ^(iLpL+L+l) ^ _ _ _ ^ f(iL P L + (Lp-l)L+0^ 



0«X J) = ' (Mi)^',f2) 

which ensures that 7-* is of full column rank. Similar to the analysis that was done for DS-OFDM systems, it is 
necessary to have 

' ' (46) 



(44) 



(45) 



L P LN T < N P 



The FDKD pilot pattern has only one non-zero pilot symbol at the center of the pilot cluster, i.e., the only 
non-zero column in Pp2 t is the (Lp — l)/2st column. If we design Pl*2t such that 

Pgth-^T" (47, 
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(48) 



it follows from the definitions of Ri fe '^ and that 

c (e^) = c (r ; (0 - 1) ) =••• = £ (r,^- 1 ^) = c (f 

which ensures that "P is of full column rank. Similar to the analysis that was done for DS-OFDM systems, it is 
necessary to have 



LN T < N P ■ (49) 



V. Miscellaneous Conditions for DS-OFDM Systems and consensus with prior art 

We now derive some additional conditions on the pilot cluster length, BEM order, number of pilot clusters that 
need to be satisfied for the full column-rank condition. Additionally, these are found to be in agreement with the 
ones derived in |6| and Q- 

It follows from ([35]! that 

Note that R| f , . . . , ~R.j Np *\ all have Lp columns. However, if the FDKD pilot pattern is used, there is only one 
non-zero column (the FDKD pilot pattern has only one non-zero pilot symbol per pilot cluster) and consequently, 
we have 



r[Rf p - l! )=l. 



r(0O = r(R[ O) ) =• 

It therefore follows (depending on the pilot-pattern design) that 

£p>r(0,)>l- (50) 

Consider d20l . It follows that r (§i) = (2B C + l)r (0;). Consequently, 

L P (2B C + 1) >r(6t) > (2B c + l)- (51) 

It follows from (|23]) and (|25l) that r ((pi) = Q and r (E) = Q. Since r (AB) < min{r (A) , r (B)}, it follows from 
(|20l that r (§i) > Q. From d5l]l, it therefore follows that 

r (§i) (2B C + 1) > Q- (52) 

We now show how some of the results derived in this paper agree with those in (6) and Q- 

• Consensus with ^ : The length of the observation cluster used in (6) is O = 2B C + 1 = 1, or in other words, 
B c = 0. From (|52]i it follows that 



L P >Q ■ (53) 



This will be called the pilot-cluster BEM condition (PC-BEM). Note that PC-BEM is similar to Lemma 1 
derived in [6]. Though this is derived for a DS-OFDM system, a similar derivation can be done for the 
DS-MIMO-OFDM system. 
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TABLE I 

Parameter Sets used in Simulation Studies. 



Parameter 
Set 


Parameter Set Values 


Si 


{A = 128, Q = 3,P b = 1, L P = 3, L? c = 1, P scp = 8, L = 4, / D = 0.1, L P L = 12, JV> = 16, QL = 12} 


s 2 


{N = 256, Q = 3,Pb = 1, L F = 3, B c = 1, P scp = 16, L = 4, f D = 0.1, L F L = 12, A P = 16, QL = 12} 


s 3 


{A = 512, Q = 3, P{, = 1, L P = 3, P c = 1, P scp = 16, L = 4, / D = 0.1, L P L = 12, A P = 32, QL = 12} 


s 4 


{A = 1024, Q = 5, P b = 2, L P = 5, B c = 2, P sop = 16, L = 4, f D = 0.3, L P L = 20, N P = 64, QL = 20} 


TABLE II 

Summary of the Scenarios for which the Full Column-Rank Conditions are tested. 


Parameter 
Set 


Remarks 


Si 


Satisfies LpL < JVp (37). Pilot patterns designed as per (34) and Sec. 1V.B. ensure lull column-rank ot P in DS-OEDM and 
FDKD-based DS-OFDM systems, respectively. 


s 2 


Satisfies LpL < Ap (37). Pilot patterns designed as per (34) and Sec. 1V.B. ensure lull column-rank ot P in DS-OEDM and 
FDKD-based DS-OFDM systems, respectively. 


s 3 


• Satisfies LpL < Ap (37). Pilot patterns designed as per (34t and Sec. IV.B. ensure full column-rank of P in DS-OFDM 
and FDKD-based DS-OFDM systems, respectively. 

• Satisfies L p LNt < Ap f46"V Pilot patterns designed as per (44) ensure full column-rank of P in a DS-MIMO-OFDM 
system with At = 2 transmitters. 

• Satisfies LNt < Ap (49t. Pilot patterns designed as per i47\ ensure full column-rank of P in a FDKD-based DS- 
MIMO-OFDM system with A" T = 2, . . . , 6 transmitters. 


s 4 


• Satisfies LpL < Ap 437 1. Pilot patterns designed as per (34) and Sec. IV.B. ensure full column-rank of P in DS-OFDM 
and FDKD-based DS-OFDM systems, respectively. 

• Satisfies L p LNt < Ap (46\. Pilot patterns designed as per (44) ensure full column-rank of P in a DS-MIMO-OFDM 
system with At = 2,3 transmitters. 

• Satisfies LNt < N p (49). Pilot patterns designed as per (47) ensure full column-rank of P in a FDKD-based DS- 
MIMO-OFDM system with At = 2, . . . , 6 transmitters. 



Consensus with ^ : Assumption 1 in states that all subcarriers are used for channel estimation (all 
subcarriers are pilot subcarriers). Assumption 2 in Q states that the pilot cluster is a FDKD pilot cluster. It 
therefore follows that N = NpLp,Lp — 2B C + 1,G = 0. Since FDKD pilot patterns are used, it follows 
from (f39b that Np > L. This satisfies a part of the condition stated in [7, (18)]. 

It follows from (O that r{B{] = 1. From the definition of f? ; in d20t . it follows that r (f? ; ) = 2B C + 1 = L P . 
Consequently, ( |52] | translates to hp > Q. Since N = NpLp, we have N/Np > Q. Together with the previous 
condition, we now have 

N 

(54) 



N 
Q 



> N P > L 



which completely agrees with [7, (18)]. 
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VI. Simulation Studies 

Four BEM models are considered in the simulation studies, namely, P-BEM, CE-BEM, S-BEM and GCE- 
BEM. Full column-rank condition of the channel estimation matrix, V, is tested in various scenarios. Though not 
explicitly mentioned, it will be understood that all the four BEM models will be tested in each scenario. All the 
four BEM models satisfy the BEMC for each of the parameter sets in Table U The RNC-BEM condition is tested 
for the parameter sets in Table [Q in the context of DS-OFDM, FDKD-based DS-OFDM, DS-MIMO-OFDM and 
FDKD-based DS-MIMO-OFDM systems. In each case, it is found that all the BEM models satisfy the RNC-BEM 
condition. 

The full column-rank condition of the channel estimation matrix, "P, is tested for the parameter sets in Table U in 

(i) 

the context of a DS-OFDM system and FDKD-based DS-OFDM system. When the pilot-pattern matrix, P pE [t, is 
designed according to d34l l and Sec. III.B, it ensures the full column-rank condition of V in DS-OFDM and FDKD- 
based DS-OFDM systems, respectively. Similarly, for the case of DS-MIMO-OFDM system, Parameter set S3 is 
tested with two transmitters (Nt — 2), while Parameter set S4 is tested with two and three transmitters (Nt — 2, 3). 
It is found that the pilot-pattern matrix, Ppl t , designed as per (l44t . ensures the full column-rank condition of "P 
in a DS-MIMO-OFDM system. Furthermore, for the case of FDKD-based DS-MIMO-OFDM system, Parameter 
sets S3, S4 are tested with two to six transmitters (Nt = 2, ... ,6). It is found that the pilot-pattern matrix, P p a t , 
designed as per (l47b . ensures the full column-rank condition of V in a FDKD-based DS-MIMO-OFDM system. 
The summary of the various scenarios are given in Table [HI 

VII. Conclusions 

We derived and discussed the conditions for the full column-rank condition of the channel estimation matrix in 
a DS-OFDM system. Firstly, it was shown that a BEM model needs to satisfy a BEMC. Secondly, the BEM model 
and the pilot-pattern design need to satisfy the RNC-BEM condition. The RNC-BEM condition was derived in the 
context of a DS-OFDM system. Once the BEM and RNC-BEM conditions were satisfied, it was shown how to 
design the pilot patterns (the pilot-pattern matrix) that ensure the channel estimation matrix is of full column rank 
in a DS-OFDM system. As special cases, we considered the FDKD pilot pattern. The pilot-pattern design method 
was extended to the case of a DS-MIMO-OFDM system with as many as six transmitters. The PC-BEM condition, 
which relates the pilot cluster length to the BEM order, was derived for a DS-MIMO-OFDM system and was shown 
to be similar to that published in a recent journal. The relationships between the number of subcarriers, number of 
pilot clusters, BEM order and multipath length was derived and was shown to be similar to that derived in a recent 
conference paper. 
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VIII. Appendix RNC-BEM 

The RNC-BEM condition will be derived only for DS-OFDM systems. For other cases like DS-MIMO-OFDM 
systems, one can derive it in a similar way. The purpose of the derivation is to show that the example of pilot-pattern 
design and BEM models considered in this paper indeed satisfy the RNC-BEM condition. Hence, we mathematically 
derive the relationship between the pilot patterns and BEM matrix for the sake of completeness though one can 
easily verify the RNC-BEM condition in Matlab or a computer-aided approach. 

Let eo, . . . , eQ_i, denote the basis vectors of E. The Q orthonormal bases that span the column space of B are 
denoted by b , . . . , bg_!. Furthermore, let Q 

w M=wWb„« = ) ... ) Q-l. (55) 
From d20b and the rank-nullity theorem, we have the relationship 

(56) 

The RNC-BEM condition states that none of the vectors eo, . . . , eQ_i, should lie in the null space of 6\, i.e., 

§ie Q ^O,q = 0,...,Q-l, l = Q,...,L-l- 

In what follows, we derive the condition for a given I and q. From ( f2Qb and d56l l. it follows that Oieq ^ implies 
that at least one of ^w^' 9 ' , i = —B c , . . . , B c , should be a non-zero vector. Assuming that w^' 9 ) ^ 0, we analyze 
the conditions when 9i-w^' q ^ ^ 0. 

Recall from d35l l that when P pa t is designed as in d34i l. we have 



A-Bo,q)T w (Bc9)Tl T , 



C 



(Rf)=-..=£(Rf p - 1 )). (57) 

Furthermore, the various columns of r| ^ are orthogonal to each other, i.e., Rj ^Rj = L p 1l p - Define an index 
vector p 3 as p 3 = [j, j + L p ,..., (N P - l)L P + j] , j = 0, . . . , L P - 1. 
We can write 

j=L P -l 

e iw ^= 0i{:Pi> {i ' q) {Pi)- (58) 

Recall from Fig. [2] that all columns of 9i(:,pj) span the same column space. From the definitions in ( [TJl i, d35l >, 
(IBTT i and Fig. |2] it follows that 

C (*,(:, Pj (0))) = ...=£ (ei(;, Pj (N P - 1))) = C (fC^+O) , < j < L P - 1- (59) 



If b q lies in the null space of WW, then wC'"' = 0. If b q does not lie in the null space of W''l, then from the rank-nullity theorem, 
w (»,9) i s a basis vector of W( l 'B. Note that if the BEMC is satisfied there is at least one non-zero w^ ,? ',i = —B c , • ■ • , B c . 
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j = 0,...,L p -l 



-B C ,...,B C 



Fig. 4. Structure of the Matrix A^' 9 '. 



The choice of P pa t as in ( [34-b ensures that 

C(O l (:,p J )) = c(fU L+t >) (60) 

which implies that 

CiOi^p^LCie^pk^J^k (61) 

i.e., 

/ (:,p J ) H ^(:,p fe ) = O,j^fc. (62) 

This means that for Oiw 1 - 1 -^ to be a non-zero vector, at least one of the summands in X^=o 1 Pi) w ^' 9 ' ) (Pi) 
should be a non-zero vector. Let A^' 9 ^ be a matrix associated with §i and eQ. Furthermore, A^ Lq \i^j) is equal to 
1 if 6i(:, pj)yv^ 1 ^ (pj) is a non-zero vector and equal to if 0i{:, Pj) r w lyl,q ' (pj) is a null vector. This is depicted 
in Fig. |U The RNC-BEM condition can be now stated as 

A ( *' 9) ^0, l = 0,...,L-l, q = 0,...,Q-l- (63) 

Note that A^' 9 ' 7^ implies that e q does not lie in the null space of §1. We now derive A 1 -- 1 '^ (i, j) in terms of b q , 
the gth basis vector of the BEM matrix B. For this we need to examine, in detail, the quantity w^' 9 * 1 (p^). 

From d55l ), it follows that w^' 9 ) (p.,-) = WW(pj, :)h q . Consider the nth column of W^(pj,:). From the 
definition of WW in ( fT3l , it follows that 

where 7<J> 1 ) = pi(j) — Pb — i- Note that <?;(:, pj) is a symmetric (circulant) matrix and as stated before, 

£(^(:,P,))=^(f (jX+0 )- (65) 

This implies that 

0,(:,Pj)W«(p;^) ^ 0, (n)^, = 3^ + 1, 

OlO,Pj)WW(p J - > n) - 0, (n) Wp ^ jli + J (66) 
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Define = [jL + IJL + 1 + N P , ...JL + 1 + (P scp - 1)N P }. From §4$, §5$ and <|66]), it follows that 



where = Efc=o 1 h i (n (i '°(fc)) e i2,r T W ' <)A/nW ' !) W. Since it follows from (US that 0;(:, p^f^+O* ^ 0, 

A^ ,9 )(i,j) is equal (not equal) to if ^i(*>«> l >^) is equal (not equal) to 0. We therefore arrive at 
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